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Positron Annihilation on Polycyclic Aromatic 
Hydrocarbon molecules in the Interstellar Medium 
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ABSTRACT 

We examine the annihilation of positrons on polycyclic aromatic hydrocarbon (PAH) 
molecules in interstellar medium conditions. We estimate the annihilation rates of 
positrons on PAHs by a semi-empirical approach. We show that PAHs can play a 
significant role in the overall galactic positron annihilation picture and use the annihi- 
lation rates and INTEGRAL galactic emission measurements to constrain the amount 
of PAHs present in the ISM. We find an upper limit of 4.6 x 10~ 7 for the PAH abun- 
dance (by number, relative to hydrogen). 

Key words: ISM: abundances - ISM: molecules - gamma-rays: theory. 



1 INTRODUCTION: POSITRON ISM 
ANNIHILATION ISSUES AND 
INTERROGATIONS 

In Astrophysics, the origin of positrons annihilating in the 
interstellar medium (ISM) and the nature and distribution 
of PAH molecules in space represent two unsolved problems 
which do not, at least on the surface, appear to have any- 
thing to do with one another. 

The abundant presence of positrons in our galaxy, 
especially in its bulge, and the high rate of annihila- 
tion determined through its gamma-ray line signature 
at 51 1 keV have been established for a few decades 
now ( Leventhal. MacCallum. fc Stangl 19781: iPurcell et al.l 
19971: lHarris et al.lll998l ; iMilne et alj|200d : I Jean et alj|2003t 
Knodlseder et al.1 l2005t). after the initial observation of 
the radiation by Johnson. Harnden. fc Havmes 1 19721 ) . 



mained mostly a mystery, for the simple reason that no 
single astronomical population of objects is known to have 
both the spatial distribution that the annihilation radia- 
tion map presents and the capability to produce such high 
rates of positrons (for a brief rev iew of the problem, see 
iGuessoum, Jean, fc Prantzosll2006T ). 

The story of the PAH molecules and their role in 
the ISM parallels that of the positrons. A few decades 
ago, the radiation emission signatures of these molecules, 
namely lines in the Infra-Red spectrum of dark nebu- 
lae, seemed to point to their abundan t presence in the 



iJohnson, fc Havmes I (|l973l ), and lHavmes et al. Il|l975h . The 
flux and profile of the annihilation line, and thus the 
rate of steady production and annihilation of positrons 
have now been measured rather precisely by several instru- 
ments, balloon -borne and satellite-carried detectors (see , 
most recently, Harris et alJ Il99 1; iKnodlseder et al.l 120051; 
IWeidenspointner et al.ll2008h . The analysis of the data on 
the gamma-ray emission (line and continuum) resulting from 
the annihilation of positrons on gas in various ISM re- 
gions has allowed for a better unde rstanding of the pro- 
cesses and condit ions of annihilation |Churazov et al.ll2005l : 
Ijean et all 120061 '). Still, the origin of positrons has re- 



Galaxy l|Gillett. Forrest, fc Merrilll Il973l and others later) 
A large effort by researchers dDulev fc Williams! [l981 ; 
iLeeer fc Pugell984lAllamandola. Tielens. fc Barkerlll985f ) 
and others later) helped identify characteristics of these 
molecules, including their structure, energy levels, and 
charge states; this was supposed to help pinpoint the pres- 
ence and distribution of specific PAH molecules in different 
ISM reg ions. The latter goal has yet to be reached (see for 
instance iRuiterkamp et al1l2005l ; ICami et al.ll2005t ). so that 
even though it is now widely believed that roughly 10 % - 
perhaps more - of the carbon in the ISM can be found in 
PAH molecules (|Tielensi ri990), definitely identifying any of 
them has remained an elusive goal. 

The connection between positrons and PAH molecules 
in earthly laboratories was made about a decade ago 
when experimental positron groups started to measure 
positron-molecule annihilation cross sections and realized 
that these increase b y several orders o f magnitudes when 
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tnat tnese increase b y several orders o t magnitudes when 
molecules are large (|Surko et al ] Il988l: llwata et al. | Il995l; 
Iwata. Greaves, fc Surkol ll99ri lB arncs, Gilbert, fc Surkol 



2005t) . A substantial amount of experimental and theoretical 
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work has been performed on this topic, and measurements 
of cross sections for a small dozen PAH and alkane molecules 
bombarded by positrons have been published, leading to 
theoretical efforts to explain the strong "Feshbach vibra- 
tional" resonanc e displayed in such interactions (see the ex- 
tensive review of lSurko. Gribakin. fc Buckmanll2005h . A few 
researchers hinted at the possible importance of such reac- 
tions in the ISM, but no serious consideration of this issue 
has been attempted heretofore. 

The closest that astrophysicists came to addressing this 
issue was when the interactio n of positrons with dust in 
the I SM was treated in detail (jGuessoum. Jean, fc Gillardl 
2005). In that work, dust was considered in the three forms 
it is usually presented in: "big grains" , "very small grains" , 
and PAHs. And because the positron-dust-grain cross sec- 
tion was assumed to be essentially geometric, "very small 
grains" were neglected next to "big grains", which them- 
selves proved to be relevant only in the hot phases of the ISM 
and, in some special conditions, in the warm ionized phase. 
It seemed reasonable then to conclude that PAHs, being 
even smaller and much less abundant, would play a negli- 
gible role in the positron's life and death in the ISM. But 
that was not correct, since the positron-PAH cross section 
was not geometrical but rather highly resonant, especially at 
very low energies/temperatures; indeed, with cross sections 
sometimes a million times larger than those of charge ex- 
change with hydrogen atoms, abundances of PAHs of about 
10 -6 (by number, Npah/Nh) would make them roughly as 
important as other species. 

We must stress, however, that PAHs, although com- 
monly considered as the "molecular end" of the dust grain 
size distribution, behave very differently with respect to 
positrons. That is why the treatment of positron annihila- 
tion on PAHs is warranted now and is fundamentally distinct 
from all past (astrophysical) works. 

In this work, we investigate the extent to which PAH 
molecules are relevant in the positrons' annihilation in 
the ISM. We have thus assembled the relevant informa- 
tion regarding PAHs in the ISM (abundances, sizes, dis- 
tributions, charge states, etc.), and on the interactions be- 
tween positrons and PAHs. With that we estimate the rate 
of positron annihilation on PAH molecules in the various 
ISM phases and compare it to those of other processes 
jGuessoum. Jean, fc Gillardll2005l '). 

In the next section we briefly review current knowledge 
about PAH molecules in the ISM. In section 3 we present 
the main experimental information on positron annihilation 
on PAHs and then proceed to determine the positron-PAH 
reaction rate at 8000 K (i.e. in the warm media) by sev- 
eral steps. In section 4 we use that knowledge to calculate 
the rate of annihilation of positrons on PAHs in the rele- 
vant phases of the ISM, paying particular attention to the 
PAH charging effect. In section 5 we explore the observa- 
tional consequences of our calculations using the measure- 
ments performed with the spectrometre SPI on the space 
observatory INTEGRAL (INTErnational Gamma-Ray As- 
trophysics Laboratory). In section 6 we summarize our find- 
ings and present our main conclusions regarding the possible 
importance of this process in current positron astrophysical 
studies and point to the work that will be needed in the fu- 
ture in order to advance our understanding of this problem. 



2 PAHS IN ASTROPHYSICS 

Polycyclic Aromatic Hydrocarbons are organic molecules 
which consist entirely of C and H atoms and have a poly- 
cyclic structure, which makes the energy needed for break- 
ing them up much higher than average; this then helps them 
survive in space. In addition, each ring has 6 electrons from 
the carbon atoms "floating" around and contributing to the 
binding. These molecules have "aromatic" properties. The 
prototype example of a PAH is Naphthalene (CioHg), which 
is made of two rings. A simpler example would be benzene 
(CsHe), but it consists of only one ring, so it is not "poly- 
cyclic" . 

These molecules became the focus of important as- 
trophysical research when the "unidentified infrared (UIR) 
emission bands" (3-13 /im) that had been observed from 
nebu l ae since the early seventies were shown (|Leger fc Puged 
1 1984 lAllamandola. Tielens. fc Barkerl Il985l ), more than a 
decade later, to be very similar to those produced by PAHs 
when they are temporarily heated (in space by UV radi- 
ation); UIR bands then came to be referred to as AIBs 
(aromatic IR bands). Indeed, the usual dust grain popu- 
lations could not be responsible for such emissions because 
they were too big and thus could not be heated to 1000 
K or more and cooled quickly. Using observations made by 
space-borne IR telescopes, similar IR emission bands were 
later found i n a variety of objects, ranging from comets to 
galaxies (e.g. lEhrenfreund et al.ll2002l ). ISalama et all (|l999h 
concluded that "PAHs are ubiquitous throughout the gen- 
eral diffuse ISM". Such identifications allowed for an infer- 
ence of the abundance of PAHs in the ISM, ro ughly believe d 
to be 10~ 7 -10~ 6 (by number, i.e. N PA h/N h ) l|Allainlll996l '). 
making them the most abundant molecules in the ISM after 
H 2 and CO. 

PAHs are often seen as an extension of, or even the 
seeds for, the sma l l dust grain populations (lAllainl 19961 ; 
ISalama et allll999l ; iPeeters II2OO2I ; lAbergel et al.1 1200a ). In- 
deed, not only are they the largest molecules k nown in space , 
but the largest of them have been described (|Peeters II2002I ; 
iRapaciolil 12005) as aggregates of planar molecules, which 
when stacked can become the "very small grains" of dust, 
which range in size from 1 nm to 10 nm. In fact, it is pos- 
tulated that PAHs constitute the first step in the formation 
of dust grains. 

The above-mentioned emission lines are generally as- 
sociated with C-C and C-H vibrationa l modes of PAHs of 
various sizes (iDulev fc Williams! [l98ll first, and many re- 
searchers subsequently) . It is also often concluded that PAHs 
with a rather large number of carbons (greater than about 
30) are responsible for some of the AIB features, particu- 
larly those around 6.2 fj,m and 11.3 /im. Since the photodis- 
sociation threshold energy is a few tenths of eV per car- 
bon atom, it is clear that "small" molecules (with a num- 
ber of carbon atoms Nc ~ 30) would more easily be de- 
pleted in H I I regions by energet ic photons (|Allainl 1 19961 : 
iPeeterslliool lAbergel et alll2005l ); bigger ones easily sur- 
vive. It has also been postulated that UV, cosmic rays, and 
shock waves could chemically alter such molecules, particu- 
lar in warm/hot, collisionally active environments. Finally, 
PAHs can be neutral or electrically charged, depending on 
the de nsities of the photoionizing radiation and of free elec- 
trons (|Omont Hl98lil ; lDraine fc Sutinll 19871 : iLepp et al.lll988l : 
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Table 1. Measured Z e f t for PAH molecules at room temperature 
(Iwata et al. 1996; Iwata et al. 1997). 



Molecule 


Chemical Formula 




Benzene 


C 6 H 6 


18000 


Toluene 


C 7 H 8 


190000 


Naphtalene 


CioHg 


494000 


Anthracene 


C14H10 


4330000 



Bakes & Tielens 1994: Allainl 


1996j; 


Da 


rtois & d'Hendecourtl 


1997; Wcineartner & Draind 


2001a 


b; 


Abereel et al.ll2005h. 



and their IR emission spectra are then significantly differ- 
ent. 



3 CURRENT LABORATORY DATA ON 
POSITRON ANNIHILATION ON PAHS 

As mentioned in the introduction, interest in e + -PAH anni- 
hilation grew with the development of experimental tech- 
niques that allowed for the measurement of cross sec- 
tions of positrons annihilat ing with more and more com- 
plex a toms and molecul e s (llwata. Greaves, fc Surkol 19961. 



1997] ; llwata et all 120001; ISurko et all |2000|; iGilbert et all 
200ol: ISullivan et all l2002l: iBarnes. Gilbert, fc Surkol 1 200^ 
Marler et all |2004| ; [Surko. Gribakin. fc Buckmanl 12005 



Surkoll2006l ). Indeed, it soon became apparent that e + -PAH 



reactions exhibited huge cross sections, or equivalently reac- 
tion rates or Z e ff, the latter being defined as: 



A 



(1) 



where A is the annihilation rate (in cm /s), and ro is the 
classical electron radius; it turns out that Z e ff has es- 
sentially nothing to do with the number of electrons Z 
of the species (for example, Z e ff — 8 for H), and in- 
deed reaches values of 10 7 or more for PAH molecules 
(see the references given above). However, only rela- 
tively small PAH molecules have been experim e ntally 
investigated so far jh yata, Greaves, fc Surkol Il996l , 1 19971 ; 
IBarnes. Gilbert, fc Surkol 120031 ). as molecules with more 
than 3 aromatic rings have such low vapor pressures that 
performing those measurements is a great experimental chal- 
lenge. The only PAH s for which & Z e ff was actually exper - 



1997) 



1996) 



imentally measured (llwata. Greaves, fc Surkcl ( l996 
are listed in Table 1. In the llwata. Greaves, fc Surkc 
experiment the FWHM of the annihilation line resulting 
from positrons on Naphthalene was also measured and found 
to be 2.29 keV (compared to 1.71 keV for H 2 ). 

The very large values for Z £ ff have been theoretically 
interpreted in terms of resonant interaction between the 
positron and the molecule, or equivalently by t he forma- 
tion o f temporary positron-molecule bound states. iGribakinl 
(2000) has come up with a theoretical model to repro- 
duce the experimental data on the cross sectional profiles 
in terms of the positron momentum distribution and of 
the temperature of the medium; he also produced a sim- 
ple fit for the exponential increase of Z^ff in terms of 
the number of atoms in the aromatic molecules (he found 
Z a ff oc A^ 8 ' 2 , where N is the total number of atoms in the 



molecule). We note, however, that Z e ff/Z tends to saturate 
for larger N's in other types of m olecules (see, for instance, 
the d ata on alkanes in Fig. 8 of IBarnes. Gilbert, fc Surkol 
(2003)). None of the studies which investigated the varia- 
tions of Z e ff with the physical characteristics of molecules 
could satis factorily reproduce the saturation effect for large 
molecules jMurphv fc Surkolll99ll; lLaricchia fc Wilkinlll997l , 
1 19981 ; llwata et al.ll2000i r This saturation is clearer when val- 
ues of Z e ///Z are presented as a function of the n umber of 
electr ons Z of the molecules (see Iwata et al.lll995l and Fig. 
39 of lSurko. Gribakin. fc Buckmanll2005l7 and we thus pro- 
pose a different fit. We have found the relation : 

, n( ^L )=j4 (l_ e -#) ( 2 ) 

to give a satisfactory fit (see Fig. 1); the number of measure- 
ments is much greater for alkanes than for PAHs; we show 
that Zejfs for the two types of molecules converge in this 
representation. In our fitting function (Eq. 2), the best pa- 
rameter values are: A = 10.75±0.44, B = 40.1±3.8. We has- 
ten to add that this relation was, like Gribakin's, not based 
on any physical considerations, only the fact that experi- 
ments show that Z e f f /Z levels off for large molecules. We 
should add that e~-PAH cross sections, being non- resonant, 
are not of much help in this regard; like Gribakin, our only 
recourse was an empirical fit. 

We must emphasize that the experimental values of 
Z e f f were obtained at room temperature only, whereas the 
temperature we are mainly interested in is 8000 K, since 
positron annihilation in the ISM is believed to occur mostly 
in the warm (neutral and ionized) phases , where T ~ 8000 K 
jChurazov et al.ll2005l ; I Jean et al.| [2006) and PAHs arc ex- 
pected to be easily evaporated in the hot phase (where T 
~ 10 6 K). To our knowledge, o nly one experiment has mea- 
sured Z e ff as a function of T (jlwata et al.ll200ol ). but this 
was limited to very few and small molecules (methane CH4, 
ethylene C2H4, and butane C4H10) and to temperatures up 
to about 2500 K. More importantly, we note that measure- 
ment errors are largely unknown (rarely, if ever, given in 
the scant literature) ; consequently, the uncertainties we have 
given above (for A and B) are related to the fit only, and 
the few measurements that are available do not allow us to 
constrain the "saturation" of the Z e f / function (in terms of 
Z) very well. We do insist, however, on the existence of such 
a saturation, unlike what the Gribakin function implies, and 
that in turn has important consequences on the relative role 
of PAH's in the annihilation of positrons in the ISM. 

So we need to determine the value of Z e ff at 
8000 K, for which we need measurements of the 
cross section as a f unction of the positron ene rgy. 
Some experiments (see Barnes. Gilbert, fc S urko 200JB] and 



ISurko. Gribakin. fc Buckmanl I2OO5F have performed such 
measurements, so we can simply calculate Z e f / as a function 
of the positron temperature by way of the reaction rate A: 



A = (ov) 



-E/kT 



/tt (kT) 



3/2 



a(E)vdE . 



(3) 



where a is the annihilation cross section, v the positron ve- 
locity, and E its kinetic energy. 

Noting th at at low tempera tures Z e ff drops as T _1//2 
(see Figure 8 of llwata et al|[2000T ). the above approach would 
yield the right behavior for Z e ff only if the cross sections 
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1 1 

20 40 60 80 100 120 140 

Number of electrons Z 

Figure 1. Comparison of ln(Z e f f /Z) for alkanes and PAHs as functions of the number of electrons in the molecule Z (Iwata et al. 1995). 
Also shown is the best fit of the empirical function (see Eq. 1) obtained with A = 10.75 and B = 40.1. 



Table 2. Calculated Z e ff at T = 8000 K for a few (medium- 
large) Alkane (top group) and PAH molecules (bottom group). 



Molecule 


Chemical Formula 


Z eff 


Butane 


C4H10 


6.4 X 10 3 


Hexane 


C6H14 


5.5 x 10 4 


Heptane 


C 7 H 16 


1.8 x 10 5 


Octane 


CgHig 


3.5 x 10 5 


Nonane 


CgH20 


6.8 x 10 5 


n-Dodecane 


C12H26 


2.6 x 10 6 


Naphthalene 


CioHg 


4.2 x 10 5 


Anthracene 


C14H10 


4.1 x 10 6 


Hexacene 


C26H16 


2.4 x 10 7 


Octacene 


C34H20 


1.1 x 10 8 


Decacene 


C42H24 


1.5 x 10 s 



were extended at very low energies as E' 1 . We were thus 
able to reproduce the experimental data of Z e ff(T) for bu- 
tane rather well for most temperature ranges, although Z e f / 
tends to decline very slowly at higher temperatures and not 
flatten out completely. Figures 2 (a and b) show the behavior 
of Z e ff as a function of temperature for various molecules; 
Figure 2b shows the variation of normalized values of Z e ff 
on a linear scale so as to emphasize the weak dependence of 
Zeff on T. 

Having noted the variation of Z e ff with temperature 
for alkanes, which we assume to be the same for PAHs, we 
then produce values of Z e f / at 8000 K for a few PAHs using 
Figure 1 and the T-dependence. In Table 2 we give values 
of Zeff for some alkane and PAH molecules. 

We emphasize the large uncertainties inherent in these 
calculations due to the fact that only a few actual measure- 
ments exist for positron annihilation reactions with alkanes 
and/or small PAHs at low temperatures. We call for new 
experiments and extensive measurements as well as theoret- 



ical calculations of annihilation cross sections for positrons 
with as many PAH species as possib le. 

iBarnes. Gilbert, fc Surkol (2003) measured the rate of 
annihilation of low energy positrons with many hydrocarbon 
molecules. They found that Z e f j is enhanced at positron en- 
ergies that correspond to resonance energies of the strongest 
infrared active vibrational modes. They interpreted these 
results as excitations of vibrational Feshbach resonances by 
incoming positrons which are temporarily trapped on the 
molecules. Based on this observation, it would in principle 
be possible to infer the variation of the annihilation rate 
on PAHs as a function of the temperature of thermalized 
positrons by integrating the product of their (Maxwellian) 
velocity distribution with the spectrum of the vibrational 
modes of the PAH molecules. However, the observed reso- 
nance effect does not seem to be systematic. Indeed, there is 
no clear enhancement of Z e ff for some energies of positron s 
(e.g. see Figs. 12 and 13 of IBarnes. Gilbert, fc Surkoll2003f l. 
Moreover, a shift in the position of the resonance with re- 
spect to the vibrational mode is observed, which seems to be 
linked to the binding energy of the positron to the molecule. 

Referring to recent work on energy-resolved 

positron-molecule ann ihilat ion cross sections 

ijBarnes. Young, fc Surkol 120061 ). Surko has remarked 
(private communication) that the temperature dependence 
of Z e ff referred to above (Zeff oc T -1//2 obtained by 
llwata et all l200d up to 2500 K) may not hold at 8000 
K, since the 1/E variation of the cross section, which 
leads to the T -1 ^ 2 dependence of Z e ff, seems to apply 
only at very low energies (E « E r es), that is for direct 
annihilation, whereas for "high" temperatures, resonant 
effects change that behavior. Should Zeff decrease faster 
with T, the values at 8000 K would be significantly smaller 
than inferred here. However, this uncertain behavior is 
itself unclear, as we do not know the binding energies for 
the larger PAHs and thus the energy difference between 
those and the vibrational energy values. This shows that 
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Figure 2. Z e tt as a function of temperature for various molecules (alkanes and PAHs) in (a) logarithmic scale and (b) normalized to 
the value 1. 



additional (unknown) factors need to be taken into account 
in order to build a correct model. 

Consequently, we cannot yet analytically obtain Z e ff 
as a function of the temperature and of N (the number of 
atoms) in large PAHs with our current level of understand- 
ing of the annihilation of positrons with molecules. We are 
then left with the necessity of producing a numerical fit of 
the available data for Z e //(N) and then extrapolating for 
larger molecules, since it turns out that the prime contribu- 
tors for positron annihilation on PAHs in the ISM are those 
with a number of atoms around 50. (See below the discus- 
sion about PAH size distributions and the contribution of 
each species to the annihilation rate.) 

Since the larger PAH molecules will be the more rel- 
evant ones for our astrophysical positron considerations, 
we believe that this issue of Z e ff for PAHs with Nc XL 30 



(N 50) is an important one that should be attempted ex- 
perimentally at the earliest. 



4 ANNIHILATION OF GALACTIC 
POSITRONS ON PAHS (RELATIVE 
GALACTIC IMPORTANCE) 

We now set out to determine the rate of positron annihi- 
lation on PAH molecules in the ISM and compare it with 
those of positrons with the usual gas components (free elec- 
trons, atomic and molecular hydrogen, helium, dust grains) . 
We may start from the reaction rate formula A = Z e f$-Kr\c 
defining Z e ff, and sum up the contributions of the various 
PAH species through their number-abundances y s (the sub- 
script s referring to any given species) , considering that some 
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may be partially destroyed, depending on the physical con- 
ditions in the medium. One must also take into account the 
possible electric charge of the molecule, which would then 
enhance or reduce the "affinity" between the positrons and 
the molecules. 

The reaction rate for each PAH species can then be 
written as: 

A e +_ s = y s Z efs , s nrlcf e i ec - s , (4) 

and this would then be summed over for the various PAH 
species. Note that y 3 includes the destruction effects; that is, 
2/s = t/so x fdest-s, where y 3 o represents the abundance of an 
s-species of PAH in the absence of any destruction effects; 
fdest-s and feiec-s will be estimated in the next sub-section. 



reaction rate gets multiplied by an "electric charge fac- 
tor", f e i e c, where f e ; ec is given by (1 - Ze 2 /a s kT) or 
exp(— Ze 2 /a s kT) depending on whether the molecule, with 
a radius a s , is negatively or positively charged, and where 
feiec then takes the values of 3.1 and 0.12 in the Z = -1/+1 
charge state, for a typical PAH size of 1 nm. 

As to the PAH destruction effect, we simply consider 
that in the hot phase of the ISM, electrons are more than 
energetic enough (kT ~ 100 eV) to break the molecule every 
time a collision occurs, so that the PAHs are largely depleted 
from such an environment; however, in the warm and cold 
phases (T fa 8000 K and 10-100 K respectively), collisions 
will rarely, if ever, destroy the PAH molecules; therefore fdest 
= for the hot phase and fdest = 1 for the others. 



4.1 PAH charging and destruction: 

Several studies have been performed on the determination 
of the charge states of dust grains and PAH molecules of 
various sizes in different cond i tions and by several physi- 
cal processes. iDraine fc Sutml (|l987T l considered collisional 
effects and showed that charge states other than 0, +1, 
and -1 have extremely low probabilities (~ 10 -4 or less) in 
various p hysic al conditions. This confirmed statements by 
lOmont I [jl986h . iLepp et all ij 19881 ) considered diffuse ISM 
clouds and determined charge states in various regions of 
the cloud (edge, center, etc.); they concluded that th e neu- 
tral state is most prevalent. iBakes fc Tielenj (|l994h then 
added photoelectric effects to the analysis and calculated 
the f(Z) probability of finding a grain/PAH molecule in a 
charge state Z. Using typical densities, temperatures, and 
UV field intensities in various phases of the ISM, they found 
(their figure 6): charge states of mostly and 1 (average 
of 0.5) and -1 and (average of ~ —0.25) for the Warm 
and Cold ISM media, respectively. The work of the se au- 
thors was revisited bv lDartois fc d'Hendeco urt (1997), who 
applied their ionization model to various conditions of the 
ISM and of nebulas; they found that in the diffuse ISM, 
65 % of large PAH molecules are negatively charged and 
35 % are neutral (different results are obtained for different 
physical conditions). P h otoelec tric charging was revisited by 
IWeineartner fc Drainei (|2001al ) with improved atomic data, 
grain size distributions, etc.; they found (their figure 9) that 
except for very large PAH molecules in the warm medium, 
where the average charge was ~ 0.4, all PAH's, especially 
mid-size molecules, will be neutral in practically all condi- 
tions. 

As to how the charge state of a PAH molecule will affect 
the reaction rate of capture and annihilation of a positron, 
we refer to the thorough discu ssion of t h e phy sics and chem- 
istry of interstellar PAHs by lOmont I ljl98rJ l. This author 
tells us that because the charge of a PAH is essentially dis- 
tributed over its entire surface, the processes which govern 
the charges of interstellar PAH's are very similar to those 
involved in the charging of interstellar dust grains. He thus 
concludes that one may reasonably adapt the classical dis- 
cussions of the charging of interstellar grains to PAH's, bas- 
ing this conclusion on several works (see references therein) . 

The positron-charged-PAH effect on the annihilation 
rate can then be dealt with in th e same way as we did in 
iGuessoum. Jean, fc Gillardl (120051). i.e. a s did lZurekl (|l985l ) 
and Guessoum. Ramatv. fc Lingenfelterl l|l99ll ). That is, the 



4.2 e+-PAH reaction rate in the ISM: 

Considering that the PAH number-abundances of 10~ 7 — 
10~ 6 quoted above are for all PAH molecules in the ISM, 
and noting (from Figure 1) that only large molecules will 
be effective annihilators of positrons, we must convolve 
the Zeff values shown in Fi gure 1 with the PAH size 
or ato m-n umber distribution; iDesert. Boulanger. fc Pugetl 
(|l99d) and [5 raine fc La zarian (Il998h give somewh at differ- 
ent PAH distributions, which iPilleri et al. I l|2009l ) present 
(their Figure 5) as a function of the number of carbon 
atoms in the molecule. The rates obtained using the two 
distributions differ by only 0.9%; in the rest of this work, 
we adopt the mo r e rec ent, "lognormal" distribution of 
IDraine fc Lazarianl (|l998l ). 

Putting all factors together then gives the following ex- 
pression for the reaction rate between positrons and PAHs: 



A 



e+-PAH 



Z eff (N)dy s (N) f ele 



■Kcr Ypah , (5) 



where Ypah refers to the total PAH abundance (relative to 
H) in the ISM. The summation over species then gives: 



A e +_ 



(6) 



We should note, however, that taking into account the 
uncertainty on the semi-empirical model of Z e f f (see section 
3) the relative uncertainty on the derived rate is about 47%. 

We now compare this result with the rates for the main 
positron processes in the various ISM phases (charge ex- 
change with H, direct annihilation with free and bound 
electrons, radiative com bination with electrons, cap t ure b y 
dust grains) as given in lGuessoum. Jean, fc Gilla rd (2005); 
the values for all rates are given in Table 3. We first note 
that if the total PAH abundance in the ISM is of the or- 
der of 10 -7 or less, then positrons will not be affected by 
PAH molecules, however large they may be. If, however, 
Ypah ~ 10 -6 , then the rate of annihilation of positrons 
by PAHs becomes non-negligible, especially if the molecules 
are neutral or negatively charged. Indeed, in the warm neu- 
tral phase this process becomes second in importance, with 
a rate up to 26 % that of charge exchange with hydrogen 
when PAHs are negatively charged. Similarly, in the warm 
ionized phase the process is then second in importance, with 
a rate up to 39 % that of radiation combination with free 
electrons. We must stress, however, that it is generally be- 
lieved that PAHs would normally exist in multiple charge 
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Table 3. Reaction rates for positron annihilation (by various processes) in the warm neutral (WNM) and ionized (WIM) phases of the 
ISM. The reaction rate for positron annihilation with PAHs is calculated using the total PAH number-abundance Ypah = 10 — 6 



Process f e +_ s (cm 3 /s) r e +_ 3 (cm 3 /s) 

WNM WIM 



Charge exchange with H 


1.8 


X 


10" 


-12 










Direct annihilation with free electrons 










1.7 


X 


io- 


-13 


Direct annihilation with bound electrons 


4.4 


X 


10" 


-14 










Radiative combination with electrons 










1.2 


X 


10" 


12 


Capture by dust grains 


6.5 


X 


io- 


-15 


4.6 


X 


iir 


-14 


Annihilation with PAH molecules 


1.5 


X 


10" 


-13 


1.5 


X 


10" 


-13 



states, with perhaps a predominance of neutral molecules - 
thus reducing their importance as annihilators. 



5 OBSERVATIONAL CONSEQUENCES 

One observational consequence of our calculation of the rate 
of positron annihilation on PAHs is the possible enhanced 
effect these molecules would have on the 511 keV line pro- 
duced in a nebula where these molecules may be abundant, 
especially if physical conditions (e.g. high collision rates) 
lead to a significant negative charging of the PAHs. And 
should PAH molecules, in such situations, start to domi- 
nate in annihilating positrons, the observed 511 keV spectra 
would be different from the usual "general ISM" ones. In- 
deed, since measurements of the s pectra of the annihilation 
of po sitrons with PAH molecules l|lwata. Greaves, fc Surkol 
1 19971 ) have determined FWHM's of the annihilation line that 
are between 2.0 and 3.0 keV, the spectra from PAH-rich 
and negatively charged nebulae would be substantially wider 
than one recorded from a warm region of the ISM, where a 
line of width 1.5 keV is produced. Furthermore, the positro- 
nium fraction (or 37/27 ratio) would likely be different in 
the two measurements. 

Using the me thod presented in 

iGuessoum. Jean, fc Gillardl (|2005l ) we have calculated 
the spectral distribution of the annihilation emission 
including the contribution of PAHs with an abundance of 
Ypah ~ 10 -6 , taking their electric charge into account. 
The FWHM of the annihilation line in PAH was here 
taken to be « 2.5 keV. In the warm neutral medium, the 
positronium fraction does not change significantly: from 
99.9 % in the absence of PAHs to 99.4 % and 98.6 % with 
neutral and negatively charged PAHs, respectively. In the 
warm neutral medium, the shape of the 511 keV line is not 
notably modified. The situation is quite different in the 
warm ionized medium. The positronium fraction decreases 
from 87 % without PAHs to 78 % and 64 % with neutral 
and negatively charged PAHs, respectively. The influence of 
PAHs distinctly deforms the base of the spectral line (see 
Fig. 4). 

Such a test cannot be undertaken at present because 
INTEGRAL-SPI does not have pinpointing capabilities to 
measure radiation from nebulas with high sensitivity, but 
the next-generation gamma-ray detectors (such as the Ad- 
vanced Compton Telescopes) should be able to perform such 
measurements. 



Another effect we have investigated is the overall con- 
tribution of PAHs to the 511 keV line from the ISM, recall- 
ing that we now have high-quality spectra of the annihila- 
tion emission fro m the galactic ce nter region measured by 
INTEGRAL-SPI l|jean et al.ll2006l ). We attempted to relate 
the reaction rate we have calculated to the annihilation spec- 
tra that have been reduced from the INTEGRAL-SPI data. 
We applied the same method as in IJean et al.1 (120061 ). but 
instead of fitting the contribution of grains (x gr in Eq. 3 of 
IJean et alfeoOrj ) we fit the parameter xpah = YpAn{feiec ) 
to the measured spectrum. We obtained an upper-limit 
xpah < 3.0 x 10 -7 . This upper-limit is in the range of 
expected values of Ypah and {feiec )■ However, the anni- 
hilation rate per unit density of PAH used for this fit is 
affected by an uncertainty of 47%. Therefore, the formal 
upper-limit should be xpah < 4.6 X 10~ 7 . This defines a 
limit in the Ypah - {feiec) space beyond which these pa- 
rameters are not allowed, or else the effect would have been 
seen in the spectral analysis of the annihilation emission 
measured by SPI. Figure 4 shows the authorized domain 
xpah = Ypah x {feiec) < 4.6 X 10 -7 . For instance, if all 
PAHs are negatively charged, then their abundance in the 
Galactic bulge should be less than 1.3 x 10 -7 . When all 
PAHs are neutral, their abundance should be less than 4.6 
x 10~ 7 . 



This result is very consistent with other estimations, 
which are obtained from totally different approaches: us- 
ing analysis of Infrared Space Observatory (ISO) observa- 
tions of the PAH emission in the Galaxy, IWolfire et al. I 
(2003) obtain a total PAH abundance of 6 x 10" 7 ; from 
their s pectral analysi s of [S iPAH]" 1 " complexes (the 6.2 am 
AIB), iJoalland etail (|2009l) infer a PAH-to-H abundance 
of 8 x 10~ 7 . We should note, however, that such measure- 
ments are not from the galactic bulge, where the presence 
and amounts of PAHs is unknown, so we cannot in all rigor 
compare the two results, despite their tantalizing closeness. 
Indeed, detections of PAH IR emission has been done from 
almost everywhere in th e G alaxy (see earlier r eferences plus 
iGiard et alll 19881 . [l989l . and lGiard et al.ll 19941 ). but not the 
bulge itself {I < 8deg); in the bulge, however, emission fea- 
tures of PAH molecules have been detected in planetary neb- 
ulae (|Perea-Calder6n et al ] 120091 : iPhillips fc Ramos-Lariosl 
120091) . 
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Figure 3. Spectral distribution of the annihilation emission in the warm ionized medium with and without PAH contribution, calculated 
assuming an abundance of PAHs of Ypah ~ 10 — 6 . The cases of neutral and negatively charged PAHs are presented separately. 
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Figure 4. Domain of Yp ah and (f e i ec ) values authorized by the SPI data: xpah < 4.6 X 10 7 . The dashed and dotted lines correspond 
to xpah = 3.0 X 10 — 7 and xpah = 1-9 x 10 — 7 , respectively. 



6 SUMMARY AND CONCLUSIONS 

The contribution of PAHs to the annihilation of positrons 
in the ISM had not been considered up till now. The simple 
reason for this is that PAHs have always been considered as 
the smallest type of dust grains, and since the cross section 
for p ositron annihilation on grains had been shown i|Zurekl 
1985) to be geometrical (modulo electric and destruction ef- 
fects), the contribution of the (extremely small) PAHs wa s 
dismissed as negligible l|Guessoum. Jean, fc Gillard! [iooBh . 
The error in that consideration was in overlooking the huge 
resonances in positron-PAH reactions (which, furthermore, 



did not occur for electrons), resonances that were vibrational 
in nature and led to increases in the annihilation cross sec- 
tions by factors of up to 10 7 . This is what prompted us to 
investigate this effect in its potential astrophysical applica- 
tions. 

Our current knowledge of e + -PAH reactions is still 
very limited, both from experimental and theoretical per- 
spectives. We have thus undertaken to estimate the rate of 
positron annihilation on PAH molecules in the ISM by a 
semi-empirical method. We have used all the information 
available to us from laboratory measurements of cross sec- 
tions and from space IR analyses and estimates of abun- 
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dances and distributions of PAHs in the ISM. We have had 
to make some extrapolations, which we have shown to be 
uncertain but rather reasonable, from alkane data to PAH 
values and from low temperatures to ISM temperatures. 

We have found that large molecules (with a total num- 
ber of atoms of about 50) are most important for positron 
annihilation. We have also obtained a total rate of positron 
annihilation on PAH molecules, and we have found that if 
the total amount of PAH (relative to hydrogen) in the ISM is 
of the order of 10 -6 , then this process becomes second in im- 
portance in the warm neutral and the warm ionized phases, 
especially if the molecules are negatively charged. The issue 
of the charge states of the PAHs in various ISM conditions 
complicates the problem somewhat; indeed, depending on 
the temperature but particularly on the intensity of the UV 
field in the region, the PAHs may be neutral, positively, or 
negatively charged, and this would enhance or reduce their 
reaction rates with the positrons. 

We have emphasized the large uncertainties inherent 
in this first attempt at tackling the problem, and we have 
pointed to the necessary future experimental, observational, 
and theoretical work for a substantial improvement of our 
understanding of this issue. In particular, we have called for 
experimental measurements of cross sections of positron an- 
nihilation on large PAH molecules. We also hope that future 
IR research, particularly with the Spitzer telescope, will help 
determine more precisely the specific PAHs which exist in 
space as well as their relative abundances in the ISM. 

Finally, we have pointed to at least one future observa- 
tional test, namely the measurement of the 511 keV line from 
specific nebulas where PAHs may be particularly abundant 
and where physical conditions may have led to substantial 
negative charging of the molecules and hence to enhanced 
positron annihilation rates. As we explained, such observa- 
tional tests cannot be undertaken with current gamma-ray 
instruments, but the next-generation gamma-ray detectors 
should be able to. We have also attempted to relate the reac- 
tion rate we have calculated to the annihilation spectra that 
have been reduced from the INTEGRAL-SPI data. We have 
shown that the spectral analysis may provide constraints on 
the abundance and the electric charge of PAH in the ISM. 
Using the semi-empirical model of annihilation of positrons 
with PAH, the analysis of the spectrum of the annihilation 
emission from the galactic bulge measured by SPI tells us 
that the number-abundance of PAH relative to hydrogen 
should be less that 4.6 x 1CT 7 if all PAHs are neutral and 
even lower if many of them are negatively charged. 
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